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Inelastic neutron scattering is used to probe the collective spin excitations of the near optimally- 
doped superconductor Ba(Fei_2,Coa:)2As2 (a:=0.065). Previous measurements on the antiferromag- 
netically ordered parents of this material show a strongly anisotropic spin-wave velocity. Here we 
measure the magnetic excitations up to 80 meV and show that a similar anisotropy persists for 
superconducting compositions. The dispersive mode measured here connects directly with the spin 
resonance previously observed in this compound. When placed on an absolute scale, our measure- 
ments show that the local- or wavevector- integrated susceptibility is larger in magnitude than that 
of the ordered parents over the energy range probed. 



I. INTRODUCTION 

The recently discovered ferropnictidc 
superconductorsi show critical temperatures ex- 
ceeding 50 K. They are interesting materials, both 
in their own right and because they provide valuable 
insights into other classes of superconductor. Theoretical 
calculations^i suggest that the electron-phonon coupling 
is too weak in the ferropnictides to produce the high 
transition temperatures which are observed experimen- 
tally. The ferropnictides have some similarities to the 
cuprates in that they are quasi two-dimensional and have 
antiferromagnetic parent compounds. Superconductivity 
in the ferropnictides can be induced^"^ from their anti- 
ferromagnetic parents by various means: electron or hole 
doping via chemical substitution; isovalent substitution 
of the iron or arsenic or by the application of pressure. 
In the widely-studied AFe2As2 (A=Ca,Sr,Ba) family 
(known as the "122" family)2r— , this can be via chemical 
substitution at the Fe or As sites or the application 
of pressure. For example, in BaFe2As2, we obtain 
superconductivity through electron doping by Co or Ni 
substitution at the Fe site; by hole doping with K at the 
Ba site; or by isovalent substitution of Fe by Ru or As 
by P. 

The proximity of superconductivity to antiferromag- 
netism suggests that the pairing mechanism in the doped 
ferropnictides is related to the spin degrees of freedom. 
Theories based on various models of the magnetic exci- 
tations have been proposed^"—. In order to take such 
theories forward, the collective magnetic excitations in 
the superconducting region of the phase diagram need to 
be characterized and the underlying interactions under- 
stood. Thus, in this paper, we report an inelastic neu- 
tron scattering (INS) study of the magnetic response, 
up to 80 meV, of a near optimally-doped supercon- 
ducting {Tc=22> K) composition Ba(Feo. 935C00. 065)2 As2. 
This composition appears not to be magnetically or- 



dered. We compare our results with similar measure- 
ments on the antiferromagnetic parent compoundsi^""— 
of the 122 series and find that (i) the magnetic excita- 
tions are stronger than in the parent compounds and (ii) 
the anisotropy in the spin-wave velocity observed in the 
parent compounds persists in this superconductor. 



II. EXPERIMENTAL DETAILS 

Single crystals of Ba(Feo.935Coo.o65)2As2 were grown 
by a self- flux method^. 15 crystals were co-aligned on a 
thin Al plate using x-rays and neutron diffraction. Our 
neutron measurements were made on an mosaic of to- 
tal mass 0.3 g. Resistivity and magnetization measure- 
ments identified the superconducting transition temper- 
ature Tc(onset) = 23 K. Elastic neutron scattering re- 
vealed no evidence of magnetic order at this doping level 
at temperatures down to 2 K. We used the MAPS in- 
strument at the ISIS spallation source. MAPS is a low- 
background direct-geometry time-of-flight chopper spec- 
trometer with position sensitive detectorsi^. A pulse of 
spallation neutrons spread over a time of about 5-10 
is produced when a pulse of protons hits a Ta target ad- 
jacent to a water moderator. Neutrons with the required 
energy (£^^=60, 80, or 140 meV in the present experi- 
ment) are then selected by an appropriately phased Fermi 
chopper (rotating at 100 Hz in the present experiment). 
The Fermi chopper is 10 m from the neutron source and 
opens for about 45 /^s. The neutrons then scatter from 
the sample (12 m from source) and are detected in posi- 
tion sensitive ^He detectors at 6 m from the sample. The 
detection time of the neutron is used to determine its en- 
ergy transfer. Data are averaged of a range of energies to 
improve the experimental statistics. The energy ranges 
of integration are given when we quote neutron energy 
transfers in the text. Data were placed on an absolute 
scale (barn sr^^ f.u.^^) by comparing the count rate with 
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FIG. 1: (color online) (a) The antiferromagnetic structure of 
the parent compound BaFe2As2. Axes refer to the tetrago- 
nal notation used throughout the paper. Solid lines denote 
the tetragonal unit cell, (b) Exchange couplings between Fe 
atoms referred to in the text, (c) Schematic illustration of 
reciprocal space (of a 2D square lattice). The ellipses repre- 
sent the anisotropic excitations reported in this paper. The 
wavevectors are used to define the spin wave velocities. 



that from a plate of vanadium. 

The magnetic cross section of an isotropic paramagnet 
is given by 
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where (7rc)^=:0.2905 barn sr~^. k; and k/ are the in- 
cident and final neutron vifavevectors and |F(Q)p is 
the isotropic magnetic form factor for a Fe^"*" orbital. 
We use Eq. [T] to convert the measured cross section to 
the energy- and momentum-dependent susceptibilityi^ 
x"(q, w). Ba(Fco.935Coo.o65)2As2 has the tetragonal 
crystal structure shown in Fig. [TJa) with lattice param- 
eters a=3.955 A and c=12.95 A. Wc use the reciprocal 
space notation to label wavevectors Q = ha* + kh* + Ic*. 
In this notation, the antiferromagnetic ordering wavevec- 
tor of BaFe2As2 is (1/2,1/2,1). The data reported in this 
paper were collected with c"^ parallel to kj. Under these 
conditions, there is a coupling of I and the u>. Wc give 
the I value corresponding to each ui in the figures and 
captions. 



III. RESULTS 

Figure [2] shows typical q-dcpcndcnt images of 
S{q,uj) = {ki/kf)(P(j/dildE in the {h,k) plane for vari- 
ous energy transfers obtained from the MAPS spectrom- 



eter. At first sight, the data look similar to those ob- 
tained over the same energy range on the related anti- 
ferromagnetic parent compound CaFe2As 2^^'"'^^ . At low 
energies, E=9.5 ± 1.5 meV [Fig. [2Ja)], the magnetic re- 
sponse is strongest at the (1/2,1/2) position. The ex- 
citations disperse with increasing energy with the pat- 
tern broadening most rapidly along the (1/2 — ^,1/24-^) 
rather than the (1/2 -|- ^, 1/2 -f ^) direction, leading to 
an cUiptically shaped response as shown schematically in 
Fig. [TJb). Fig. [3] shows cuts through the data shown in 
Fig. [2J again illustrating the anisotropic dispersion and 
broadening of the response with increasing energy trans- 
fer. 

Although Ba(Feo.935Coo.o65)2As2 is not magnetically 
ordered, its response is reminiscent of the magnetically 
ordered 122 parent compounds. In order to make a more 
quantitative analysis of our data we fitted a phenomeno- 
logical spin- wave cross section with a damped-harmonic- 
oscillator (DHO) line shape in energy. The DHO cross- 
section is specified by Eqs. 
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In these equations, W|| and vj_ describe the dispersion 
along (1/2 + ^, 1/2 + ^-type and (1/2 - ^, 1/2 + ^-type 
directions respectively (or gy and q± in Fig. [T]), F is the 
damping parameter and A is a gap (see later). Eqs.[2H3] 
give a good description of the data at all energies. The 
q-averaged or local susceptibility x"(w) is defined as: 



x"(w) 
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where the average is over sufficient Brillouin zones in the 
extended zone scheme to sample representative wavevec- 
tors of x"(q, w). The local quantity gives an indication 
of the overall strength of the excitations and can be used 
to compute the corresponding fluctuating moment 
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Figs. |3][5] show the dispersion and the energy dependence 
of determined from fits to the data. The different 

symbols in Figs. |4] and [5] represent different I values for 
the c-axis momentum. In the present experiment we av- 
erage x"(Q: ^) over < h < 1 and < fc < 1 and / over 
approximately ±0.25 to compute the values in Fig.[SJ All 
points appear to follow the same trend and thus we are 
unable to discern any evidence for magnetic coupling in 
the c-direction for the energy scale of the present exper- 
iment £^ > 10 meV. 

A second function which has been used to describe the 
magnetic response in nested Fermi liquids is the modi- 
fied Lorentzian or Sato-Maki (SM) function^i. This func- 
tion has been successfully used to model chromium2ii2^, 
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FIG. 2: (color online) (a)-(e) Constant-energy slices through the magnetic excitations in Ba(Feo.935Coo.o65)2As2 as observed 
on MAPS, (f)-(j) Fits to a phenomenological spin-wave cross-section (Eqs. [2J3]). (k)-(o) Fits to a phenomenological Sato-Maki 
cross sections (Eq.[S]) [e.g. for (n) k. — 0.21 ± 0.09, S = 0.6 ± 0.2, A = 5 ± 2]. The incident energies used were 60 meV [(a)-(b)] 
and 140 meV [(c)-(e)] and the corresponding / values were 1=1 [(a)-(b)], 2.5 (c), 4(d) and 5.5 (e). The i5i = 140 meV data were 
collected using a proton beam current of 175 nA for 80 hours. A constant background has been subtracted from each plot. 



La2_2;Srj;Cu04^ and more recently FcTei-xSear^. The 
Sato-Maki cross-section is of the form: 



resonance 
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with 



^(q) = ^{[ih-hor + ik-kof-s'y 

+ ^[{h~hof~{k-kof] 

where Qo = {ho, fco) is the nearest reciprocal lattice point 
to Q with odd {h + k). The location of the peaks in the 
response is controlled by (5, and its shape by k and A. 
For example, to obtain a peak centered at Q=(l/2,l/2), 
6 — 1/ v^. We found that the Sato-Maki function could 
give a good description of our data except at the highest 
energies £' w 80 meV. Fits of the SM-function arc shown 
in Figs.lJEl The values of x"(w) computed from the SM 
fits are statistically indistinguishable from those obtained 
using the SW cross-section. 



IV. DISCUSSION 

The magnetic response has been studied by neu- 
tron scattering in a number of superconducting ferro- 
pnictides^i"— including Ba,{Fci-xCox)2As-2r^^. Previ- 
ous experiments on the superconducting 122 systems 
have either been on powder samples^ or performed 
over a smaller energy range^i^ than the present exper- 
iment. One of the early results from INS on the iron- 
based superconductors was the observation of a 'spin- 



Thc spin-resonance is most easily in- 
terpreted in terms of a sign difference in the supercon- 
ducting gap A(k) between different parts of the Fermi 
surface^"— and is consistent with a gap function. At 
the lowest energies probed by the present experiment we 
are able to observe the spin- resonance. Fig. [5] shows this 
most directly: There is a peak in the local susceptibil- 
ity x"('^) near 10 mcV at T = 7 K which is suppressed 
on raising the temperature to 26 K. The resonance is 
strongest near q = (1/2,1/2) as shown in Fig. [H^a) 
and Fig. [31[d,h). As the energy is increased above the 
resonance energy, magnetic excitations disperse in an 
anisotropic manner (see Fig. |4]). Thus our results are 
consistent with an upwardly dispersing mode which is 
strongest near q = (1/2,1/2). This is at least quali- 
tatively consistent with theoretical predictions for a 
stated. Interestingly, this contrasts with the behavior 
in YBa2Cu306+2:, where the resonance mode disperses 
downwards in energy^^— . 

An important result from our experiment is the ob- 
servation of dispersive anisotropic spin fluctuations up 
to 80 meV. Neutron scattering studies on the anti- 
ferromagnetic parents of the 122 series, BaFe2As 2^"^'^^ 
and CaFe2As2^'ii, have observed spin-wave excita- 
tions up to about 200 meV. Unfortunately, single 
crystal measurements up to 80 meV only exist for 
CaFe2As2. The magnetic excitations have been analyzed 
using a spin-wave model based on localized moments. 
Ba(Feo.935Coo.o65)2As2 is close to magnetic order at low 
temperature. In order to parameterize our data we use a 
damped spin-wave cross-section with an energy gap. In 
our case, the gap A in Eq.[3]is due to the superconductiv- 
ity (and the formation of the resonance described above). 
We find that Eqs. [2H3] provide a good description of the 
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FIG. 3: (Color online) Constant hu and Z cuts along (1/2 — 
ft, 1/2 + ft, Z) and (1/2 - ft, 1/2 + ft, (see Fig. [ij for vari- 
ous values of hw. The blue solid line a fit to the SW cross- 
section Eqs. [2J3] and the red dotted line to the Sato-Maki 
model (Eq. [6|. Where only one line is drawn, the fits are 
indistinguishable . 



data (as illustrated by Figs. EHS] ) with T/E = 0.15 (as 
per Ref . [iTh . This method of analysis produces spin- wave 
velocities of wy = 580 ± 60 meVA, wj. = 230 ± 30 meVA 
and A = 10 ± 0.5 meV. Using the published exchange 
constants, we obtain spin-wave velocities for the parent 
compound CaFe2As2 of W|j=513 meVAi^ or 494 meVAii 
and z;±=370 meVAi^ or 348 me V Ail. If the magnetic in- 
teractions are described using a Heisenberg Hamiltonian 



<]k> 



Sfc + J2 ^ 



(7) 



where the sums are over nearest-neighbor and next- 
nearest-neighbor pairs (see Fig. [1]), we can convert the 
spin- wave velocities to effective exchange constants based 
on spin-wave theory in the ordered parent antiferromag- 
nets. Using the relations v\\ = S^/2a{Ji + 2J2) and w_l = 

SV2a^/{2J2 + Ji)(2J2 - Ji) (see Refs. [lUinill, we 

obtain Ji = 43 ± 7 and J2 = 30 ± 3 meV. It should be 
noted that this is not the only possible analysis, but it 
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FIG. 4: (Color online) The dispersion of the magnetic exci- 
tations in Ba(Feo.935Coo.o65)2As2 for T = 7 K as determined 
from fitting the damped spin-wave cross section (Eqs. [SHl]) 
along qy (a) and q^ (b) compared to that found in CaFe2 As2ii 
(blue dashed line) in each case. Symbols denote I values: ■ 
{I even); • {I odd); ♦ {I non-integer). 



is the simplest. It may also be possible to interpret the 
data in terms of the Jia-Jib-J2 modeliii^ which has been 
used to describe the ordered parent CaFe2As2. However 
\Jia - Jib\ is certainly less for Ba(Fco.935Coo.o65)2As2. 

In addition to comparing the spin- wave velocities with 
the parent compounds, wc can also compare the strength 
of the magnetic response. Fig. [5] shows the wavevector- 
averaged or local susceptibility x"(w) defined by Eq. |4l 
In the case of Ba(Feo. 935C00. 065)2 As2, we have averaged 
over < ft- < 1 and < fc < 1, and ±0.25 around var- 
ious odd, even and non-integer I values (see caption to 
Fig. [5]) . As mentioned above, points with different I ap- 
pear to follow the same trend, suggesting that there is no 
/-dependence to this partially averaged quantity. Thus, 
the graph represents the true x"(w). For comparison, we 
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Ba(Fe-| 935COo.o65)2AS2 




Energy (meV) 

FIG. 5: (Color online) The local susceptibility x"(^) at 7 K 
(closed symbols) and 26 K (open symbols). Symbols denote 
I values: ■ {I even); • {I odd); ♦ {I non-integer). Dashed line 
shows the local susceptibility for CaFe2As2^'2£. 

have computed x"{oj) for C aFe^ As^ using the exchange 
constants and Scs from Ref. llTI . This result is shown as 
the dashed line^S in Fig. [S] We note that the response is 
slightly larger in Ba(Feo.935Coo.o65)2As2 over the energy 
range investigated here. The same phenomenon occurs 
in the cuprates where, for example, the spin fluctuations 
arc stronger in optimally-doped La2-a;Sra;Cu04 than in 
La2CuO/t2^i^. One might understand this increase at 
lower energies as being due to a shift in spectral weight 
from higher energies (above the window of the present 
experiment) and the loss of the peak due to magnetic 
order present in the ordered compounds. 

It is interesting to compare the present results 
with those obtained on the cuprates. Firstly, we 
note that the magnitude of the local susceptibil- 
ity is similar in optimally doped La2_a:Sra;CuO/t2^i^ 
and Ba(Feo. 935C00. 065)2 As2. Both the cuprates and 
Ba(Feo.g35Coo.o65)2As2 show a strong response near the 
Q=(l/2,l/2). In La2-^Sr^Cu04 and YBa2Cu306+a;, 
various dispersive modes and 'resonance' features^^i^"— 
are observed below E < 50 meV. At higher energies, 
50 < < 100 meV, the response in YBa2Cu306+a- 
and La2_a;Sra;Cu04 is quasi-isotropic or has fourfold 
symmetry about the (1/2,1/2) position^^i^"— . In con- 
trast, the response in Ba(Feo.935Coo.o65)2As2 is more 
anisotropic being broader along (1/2 — ^,1/2-1-^) than 
the (1/2 -I- ^, 1/2 -I- ^). We discuss possible origins of this 
anisotropy below. 

Discussion about the best way to describe the magnetic 
interactions in the ferropnictides continues^^^'^'^. 
These materials are not Mott insulators with localized 
spins which can be described purely by near-neighbor su- 
perexchange. Neither are the ferropnictides weakly corre- 
lated metallic systems. It appears that Hund's rule cou- 
pling is responsible for producing the atomic moments 
and there is also strong magnetoelastic coupling. Keep- 



ing this is mind, we discuss our results in terms an effec- 
tive Heisenberg coupling and a weakly correlated metallic 
picture. 

There is evidence that the ferropnictides develop an 
underlying electron nematic phase^i^i^ as they ap- 
proach the magnetic ordering temperature. That is, they 
show fluctuating magnetic stripes based on the ordering 
shown in Fig. (Ha). While Ba(Feo.935Coo.o65)2As2 does 
not appear to order magnetically, it is close to the anti- 
ferromagnetic quantum critical point^^— , therefore simi- 
lar considerations should apply. It has been known for 
some timo^ that 2D square-lattice systems described by 
the Hamiltonian in Eq.[7]develop a stripe-like order of the 
type observed in the ferropnictides for J2 / Ji > 1/2 which 
is the condition here (our measurements give J2/J1 = 
0.70±0.09). Ahhough the Hamiltonian in Eq.[7]has a C4 
lattice symmetry, we expect low-frequency correlations to 
develop at low temperature through the so-called 'order 
out of disorder' mechanisn*^2^i^. This means instanta- 
neous collinear order described by two interpenetrating 
Neel sublattices as shown in Fig.[TJa). The higher-energy 
magnetic excitations of such a state with long-lived cor- 
relations of this type might be expected to be similar 
to the excitations of the corresponding ordered antiferro- 
magnct. This could explain the similarities between the 
magnetic excitations in the 10 < < 80 meV range in 
Ba(Feo.935Coo.o65)2As2 and CaFe2As2i^iii. 

The second approach which has been used to describe 
the ferropnictides is that based on nesting between the 
electron and hole pocketsiSii^ii^i^d^. In such a pic- 
ture, the magnetic excitations are correlated electron- 
hole pairs. The canonical example of an itinerant an- 
tiferromagnet which is described by such a picture is 
elemental chromiun*^. Model Lindhard calculations 
of the wavevector-dependent susceptibility xi'i,^)^'^ 
based on the band structure appear to reproduce the 
q-anisotropy in x"(q, w) observed here at higher en- 
ergies. In particular, the response is broader along 
(1/2 - ^, 1/2 + ^) rather than the (1/2 + 1/2 + ^) 
for a given energy. It is interesting to note that highly 
structured magnetic excitations, characteristic of nest- 
ing, have recently been observed in FeTci^^Se a.^^i^'''^^ . 
For FeTco.5iSeo.49, an anisotropic response similar to the 
one reported here is observed. The authors of Ref. [2^ 
find that the Sato-Maki cross section (Eq. [S]) provides 
a good description of their data over a 10-120 meV en- 
ergy range. Motivated by this we also fitted our data to 
the Sato-Maki form (see Fig. While the Sato-Maki 
function provides a reasonable description of the data, 
the phenological spin-wave cross section provides a bet- 
ter description at higher energies. 



V. CONCLUSIONS 

In summary, we have used inelastic neutron scattering 
to probe the collective spin excitations of near optimally 
doped Ba(Fei_2;Co2;)2As2 (a;=0.065). Strongly disper- 
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sive spin fluctuations are observed up to 80 meV. In the 
superconducting state, our measurements are consistent 
with a mode dispersing from the spin resonance near 
Q = (1/2,1/2) to liigh energy. At higher energies, we 
observe excitations which are anisotropic as a function of 
in-plane wavevector. The response is centered on the M 
or (1/2,1/2) position of the Brillouin zone and is broader 
along (1/2 - ^, 1/2 + C) rather than the (1/2 + ^,1/2 + 
direction. This anisotropy may be understood in terms 
of multiple J1-J2 exchange interactions or response func- 
tions derived from the band structure. When placed on 
an absolute scale, our measurements show that the local- 
or wavevector-integrated susceptibility is generally larger 



in magnitude than that of the ordered parent antiferro- 
magnets over the energy range probed here. 

Note Added. Diallo et al.— have recently reported 
magnetic fluctuations in the paramagnetic state of 
CaFe2As2 which are similar to those presented here. 
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